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ABSTRACT 

We present here a comprehensive model for the evolution of the magnetic properties of habitable 
terrestrial planets (Earth-like planets, M ^ 1 and super-Earths M ^ 1 — 10 M©) and their effects 
on the long-term protection against the atmospheric erosive action of stellar wind. Using up-to-date 
thermal evolution models and dynamo scaling laws we predict the evolution of the planetary dipole 
moment as a function of planetary mass and rotation rate. Combining these results with models for 
the evolution of the stellar wind, stellar XUV fluxes and exosphere properties of highly irradiated 
planets, we determine the properties of the planetary magnetosphere and the expected scale height 
of the atmosphere that together determine the level of thermal and non-thermal atmospheric mass 
losses. We have used this model to evaluate the early magnetic protection of the Earth and the already 
discovered potentially habitable super-Earths GJ 667Cc, Gl 581d and HD 85512b. We confirm that 
Earth-like planets, even under the highest attainable dynamo-generated magnetic field strengths, will 
lose a significant fraction of their atmospheres or their content of critical volatiles (e.g. H2O) if they 
are tidally locked in the HZ of dM stars. Planets in this mass-range with N/O-rich atmospheres, 
even under the best conditions of magnetic protection, will probably lose their atmospheres or their 
water content if they are in habitable zones closer than ~ 0.8 AU (Af^, < 0.7— 0.9Mq). Super-Earths 
Mp > seem to have better chances of preserving their atmospheres even if they are tidally 

locked around dM stars. Under similar conditions of thermal and magnetic field evolution there seems 
to exist a planetary mass-dependent inner limit inside the HZ itself, below which large atmospheric 
mass-losses in super-Earths are expected. With the nominal value of the physical parameters in our 
conservative model this limit is, for example, ^0.1 AU for 4M0 and ^0.04 AU for 8Mq. Under these 
conditions we predict that the atmosphere of GJ 667Cc has probably already been obliterated and it 
is presently uninhabitable. On the other hand, our model predicts that the atmospheres of Gl 581d 
and HD 85512b would be well protected by dynamo-generated magnetic fields even under the worst 
expected conditions of stellar aggression. 

Subject headings: Planetary systems - Planets and satellites: atmospheres, magnetic fields, physical 
evolution - Planet-star interactions 
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1. INTRODUCTION 

The discovery of extrasolar habitable planets is one of 
the most ambitious challenges in the exoplanetary re- 
search. At the time of writing, there are almost 760 
confirmed exoplanets including 55 classified as Earth- 
like planets (EPs, M ^ 1 Mm) and sup er-Earths (SEs, 
M - 1 - 10 Me, iValencia et all 120061 ). Although the 
composition of most of these planets is unknown, many 
of them should have similar compositions to the Earth 
which would make them the first extrasolar terrestrial 
planets (TPs) discovered so far. 

Among these low mass planets already discovered 
there are three c onfirmed SEs, GJ 667Cc, Gl 581d 
and HD 85512b (ISelsis et all \m)Ti IPene et al.l [20TTI : 
iKaltenegger et al] 120 111) and tens of Kepler candidates 
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(jBorucki et al.ll2011l: iBatalha et al.ll2012D that are close 
or inside the habitable zone (HZ) of their host stars If 
we include the possibility that extra-solar-system giant 
planets could harbour habitable exomoons, the number 
of already discovered potentially habitable planetary en- 
vironments beyond the Solar S ystem co uld be ri sed to 
several tens ()Kalteneggeill20ldl : [Underwood et al.ll2003i) . 
Moreover, the case for the existence of a plethora of 
other TPs and exo moons in the Galaxy is rapidly gaining 
evidence (jBoru cki et al.ll20lTI : iCatanzarite fc Sh"aoll2011l : 
iKipping et al.ir2012l ) and the chances that a large num- 
ber of potentially habitable extrasolar bodies could be 
discovered in the near future are encouraging. 

The question of which properties a planetary environ- 
ment needs in order to allow the appearance, evolution 
and diversification of life has been extens ively studied 
(for r ecent reviews see iLammer et al.l 120091 and iKastingl 
1201 0*) . Two basic and complimentary physical conditions 
must be fulfilled: the presence of a n atmosphere and th e 
existence of surface liquid water (jKasting et al.l |1993|) . 
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However, the fulfilment of these basic conditions depends 
on many complex and diverse endogenous and exogenous 
factors (for a comprehen sive enume ration of these factor s 
see e.g. iWard fc Brown lcc 2000 or lLammer et al.ll2010t ) 
The existence and long-term stability of an intense 
planetary magnetic fi eld (PMF) is one of t hese addi- 
tional factors (see e.g. iGrieBmeier et al.ll2010l and refer- 
ences therein). It has been shown that a strong enough 
PMF would protect the atmosphere of potentially hab- 
itable planets, especially its valuable content of water 
and other volatilcs, against the erosive action of the s tel- 
lar w i nd dlammer et al. 2003|, 120071 : iKhod achcnk o etHI 
l2007t iChaufrav et al. 2007.) . Planetary magnetospheres 
would also act as shields against the potentially harm- 
ful effects that the stellar and galactic cosmic rays (CR) 
could pr oduce in the life-forms evolving on its surface 
(see e.g. iGriefimei er et al.1l2005t ). Even in the case that 
life could arise and evolve on unmagnetized planets, the 
detection of atmospheric biosignatures would be also af- 
fected by a higher flux of stellar and galactic CR, espe- 
cially if the planet is around very active M-dwarfs (dM) 
(jGrenfeU et al.ll2007t ISegura et al.ll2010( ). In summary, 
a PMF would not only protect the atmosphere of the 
planet and the life growing on its surface, but also give 
us the possibility to confirm the habitability of future 
discovered planets around dM stars. 

But find suitable conditions to have TPs with strong 
enough PMFs in the HZ of the most abundant stars 
seems more problematic than previously thought. It 
has been recently predicted that most of the T Ps in 
our Galaxy could be found around dM stars (lB _ os£ 



20061 : iMavor fc Udrvl[200l iScalo et ai][2007HRauer et al 



201 It) . Actually ~ 20% of the presently confirmed super- 



Earths belong to planetary systems around stars of this 
type, including Gl 581d one of the best candi date for 
habitability presently known (jSelsis et al.l 120071 ). Plan- 
ets inside the HZ of low mass stars ( M^, < O.GMc?, ) wou ld 
be tidally locked poshi et al.lll997t iHeller et al.ll2011[ ) a 
condition that poses seri ous lim itations to their poten- 
tial habitability (see e.g. iKite et al. 2011 and references 
therein) . Tidally locked planets inside the HZ of dMs will 
have periods in the range of 5 — 100 days, a condition has 
commonly been associated with the almost complete lac k 
of a protective magnetic field (jGriefimeier et al.ll2004[ ). 
However, the relation between rotation and PMF prop- 
erties, that is critical to assess the magnetic protection of 
slowly ro tating planets, is more com plex than previously 
thought (jZuluaga fc Cuartasll2012f ). In particular a de- 
tailed knowledge of the thermal evolution of the planet 
is required to predict not only the intensity but also the 
regime (dipolar or multipolar) of the PMF for a given 
planetary mass and rotation rate. 

Several authors have extensively studied the protection 
that intrinsic PMF would provide to extrasolar plan- 
ets (IGriefimeier et al.l 120051: iKhodachenko et all 120071: 
iLammer et al.ll2007t IGriefimeier et al.l 120091 \20m . In- 
dependently the eff'ects of X-Ray and EUV radiation 
(XUV) on the thermal escape processes in weakly magne- 
tizcd planets have also been studied (Kulikoy et al.ll2006l: 
Xammcr et al. 2007, 2009; Tian et al l 120081: iTian 2009; 
Sanz-Forcada et al..,2010 '). All these works have however 



neglected the evolving nature of the PMF and have sys- 
tematically predicted field intensities using scaling-laws 
that have been revised in the last couple of years (see 



lChristensenll2010l and references therein). More impor- 
tantly the role of rotation in determining the PMF prop- 
erties that is critical in assessi ng the case of tida lly locked 
planets has been overlooked (jZuluaga fc Cuar tas 2012). 

In this work we present a model of the evolution of the 
magnetic protection of potentially habitable TPs around 
GKM main sequence stars. To achieve this goal we in- 
tegrate in a single framework the most r ecent thermal 
evolution models for this ty pe of planets ()Gaidos et al.l 
l2010l :lTachinami et al."201n. up-to-date dynamo scaling- 
laws (jChristcnscn 2010; Zuluaga & Cuartas 2(112), mod- 
els for the evolution of the stellar wind and XUV 
luminosity of low mass- stars (jGriefimeier et all 120101 : 
iSanz-Forcada et all 1201 Ij ) and the most recent results 
describing the expansion and hydrodynamical escape 
of volatiles fr om highly irradiated atmospheres of low 
mass planets (jKuhkov et al.ll2006t iTian et al.ll200l jTia^ 
l2009f l. This is the first time that all these pieces have 
been put together to produce a global picture of the mag- 
netic protection of potentially habitable TPs. 

But this is not only a model integration effort. Sev- 
eral novel features have been added to our comprehen- 
sive model: 1) we include a new treatment of the role 
of rotation in determining the PMF properties, eespe- 
cially important in assessing the magnetic protection of 
tidally locked planets, 2) we propose a phenomenologi- 
cal formula to estimate a magnetic- constrained thermal 
mass-loss rate from atmospheres protected with a strong 
PMF, and 3) we address the magnetic protection of al- 
ready discovered habitable planets and compare it with 
the case of an early magnetized Venus and the Earth in 
its current state. 

This paper is organized as follows: section [2] describes 
the model used here to calculate the evolution of mag- 
nctosphere properties. For that purpose we use a set 
of analytical fits of recently published models for the 
thermal evolution of TPs (section 12. ip and up-to-date 
power-based dynamo scaling laws (section 12. 2p . Several 
stellar key properties are required to calculate the evo- 
lution of the magnetopshere (stellar wind properties, HZ 
distances and tidally locking limits). The way we ob- 
tain these properties and their evolution are described 
in section 12.31 Atmospheric expansion models and the 
proposed phenomenological formula to model the atmo- 
spheric mass- loss rate are described in section [3] Section 
IDpresents the results of applying this model to hypothet- 
ical and already discovered habitable TPs. A discussion 
about the hypothesis on which the model relies and the 
possible sources of uncertainties are presented in section 
ini Finally several conclusions and the future prospects 
of this work are summarized in section |6l 

2. A MODEL FOR AN EVOLVING MAGNETOSPHERE 

The interaction between the PMF, the interplanetary 
magnetic field (IMF) and the stellar wind creates a mag- 
netic cavity around the planet known as the magneto- 
sphere. Magnetospheres are very complex systems but 
its basic pro perties are continuous funct ions of two ba- 
sic variables (jSiscoe fc ChristopheHll975t) : the planetary 
magnetic dipole moment A4 and the dynamical pressure 
of the stellar wind Psv/'- 
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Psw = fnnvl^ + 2nkBT. (2) 

Here i3p'P(r) is the dipolar component of the field as 
measured at distance r from the planet center, m and n 
are the typical mass of a wind particle (mostly protons) 
and its number density, respectively. Wcg = {^sw ^^^pY^^ 
is the effective average velocity of the stellar wind as 
measured in the reference frame of the planet whose or- 
bital velocity is Vp. And T is the local temperature of 
the plasma. fiQ = An x 1Q~^ H/m and are the vacuum 
permeability and Boltzmann constant respectively. 

There are three basic properties of planetary magne- 
tospheres we are interested in: 1) the maximum magne- 
topause field intensity B^p, a proxy to the fiux of high 
energy particles into the magnetospheric cavity, 2) the 
standoff or stagnation radius, Rs, a measure of the size 
of the dayside magnetosphere, and 3) the area of the po- 
lar measures the total area of the planetary 
atmosphere exposed to open field lines through which 
particles can escape to interplanetary space. The value 
of these three quantities provides information about the 
level of exposure that a habitable planet has to the ero- 
sive effects of stellar wind and the potentially harmful 
effects of the CR. 

The maximum value of the magnetopause field inten- 
sity Snip is estimated by equating the magnetic pressure 
^inp = ^mp/(2A'o) and the dynamical stellar wind pres- 
sure Psw (eq. E]) , 

Smp = {2tio)'/^P^f (3) 

Although the magnetopause fields arise from very com- 
plex processes (Chapman-Ferraro and other currents at 
the magnetosphere boundary), in simplified models i?mp 
is parametrized as a multiple of the planetary field in- 
tensity Bp as measured at the substellar point r — Rs 
(|Meadil964 iVoigt.,1995i) . 

Pmp = 2/oBp(r = Rs) = (^^) V2MRf (4) 

where /q is a numerical enhancement factor of order 1 . 
We are assuming here that the dipolar component of the 
intrinsic field dominates at magnetopause distances even 
in slightly dipolar PMF. Combining equation |3] and |4] we 
obtain an estimate of the standoff distance: 

Rs-{^^yM'/'P^-'/' (5) 

It is important to stress here that the Rs given by 
eq. ([5]) assumes that the pressure exerted by the gasses 
trapped inside the magentosphere cavity is negligible. 
This is a good approximation only if the planetary mag- 
netic field is very intense or the stellar wind is weak or 
the planetary atmosphere is not too bloated by the XUV 
radiation. In the case when any or none of these condi- 
tion are fulfilled we will refer to the Rs estimated with 
eq. ([S]) as the magnetic standoff distance. 

The last but not least important property in which we 
are interested is the area of the polar cap. The polar 
cap is the region in the magnetosphere where open field 



lines c ould transport ions int o or from the interplanetary 
space. iSiscoe fc ChenI (|1975[ ) showed that the area of the 
polar cap Ape scales with the dipole moment and the 
dynamical pressure of the stellar wind as: 

where Ms, = 7.768 x lO^^Am^ and Psw© = 2.24 nPa 
are the present values of the dipole moment of the 
Earth and the average dynamic pressure of the solar 
wind as measured at the Earth distance (jStacevI [19921 : 
IGriefimeier erall[2005h . 

In order to model the evolution of the magnetosphere 
properties we need to calculate reliable values of the sur- 
face dipolar component of the PMF B^'^p, the average 
number density n, the velocity and the temperature 
T of the stellar wind. These quantities depend in general 
on time and also on different planetary and stellar prop- 
erties. In the following sections we describe the models 
used in this work to calculate the evolving values of these 
critical quantities. 

2.1. Planetary thermal evolution 

We assume here that the main source of a global PMFs 
in TPs is the action of a dynamo powered by convection 
in a liquid iron core. Other possible sources of PMFs, 
dynamo action in a mantle of ice, water or magma or 
stellar induced magnetic fields, are not considered here. 

The properties and evolution of a core dynamo will 
depend on the internal structure and thermal history of 
the planet. Two recent works have built and solved de- 
tailed interior and thermal evolution models intended to 
study the generatio n of intrinsic PMF in su per Earths 
(|Gaidos et all IMot iTachinami eTlII [20Tl[ ) . These 
works have paid attention to differ ent and complimen - 
tary aspects of the problem: while iGaidos et al.l (|2010( ) 
have conce ntrated on the thermodynamics of the core, 
[Xachinami et al.l ()2011l ) have developed a detailed treat- 
ment of mantle rheology and convection. From these 
works a robust picture of the thern ial evolution of SEs 
is star ting to arise. More recently iZuluaga fc CuartasI 
(|2012l) have shown how by combining the results of these 
thermal evolution models and the rotational properties 
of the planets it is possible to predict not only the mean 
intensity of the magnetic field, but also its regime (dipo- 
lar or multipolar). This last property is important in 
predicting the intensity of the dipolar component of the 
field at the planetary surface and from there the dipole 
moment of the planet. 

There are three key properties that should be predicted 
by any thermal evolution model in order to calculate the 
magnetic properties of the planet: 1) the total convec- 
tive power Qconv providing the energy to be dissipated 
through dynamo action, 2) the radius of the solid inner 
core Pic and thus the height D of the convecting shell 
where the dynamo action takes place {D = Rc — Pic) 
and 3) the total dynamo life-time idyn- Figure [T] shows 
th e values of the se quantities as predicted by the model 
bv IGaidos et al.l (|201Q) for planets with the same com- 
position as the Earth in the mass range 1.0-4.8 M^. 

We have selected the model results obtained for plan- 
ets with mobile lids and habitable surface temperatures 
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Fig. 1. — Results of the thermal evolution model bv lGaidos et al.l 
Il20ig ) for TPs with the same composition as the Earth. Upper 
panel, radius of the solid inner core. Middle panel, total convec- 
tive power. Lower panel, time of the inner core nucleation (red) 
and dynamo lifetime (blue). Numerical results are indicated as 
crosses joined by dashed lines. Continuous lines correspond to the 
analytical fit used in this work to interpolate and extrapolate the 
model results (see text). 



rsurf=288 K. We assume that the potentiaUy habitable 
planets described with our magnetic protection model 
fulfill this condition at least for the time during which 
they have a dense enough atmosphere and an available 
inventory of water. 

In this model all planets start their thermal history 
with a liquid iron core with radius Re- The iron core 
cools down through the release of heat stored during the 
early phases of the formation of the planet (sensible heat) 
and generated by the decay of radionuclides. The con- 
tribution of tidal heat in this process has, however, been 
negelected. When the temperature at the center of the 
planet has decreased below a critical level, an inner solid 
core starts to form and suddenly new sources of energy 
arise (release of latent heat and buoyancy generated by 
the release of light elements). This change produces a 
strong "rebounce" of Qconv after the initiation of inner 
core growing (see middle panel in figure [T]). However, 
this phenomenon is restricted to planets with masses 
Mp < Merit » 2 - 2.5Mffl (a thresh old predicted in- 
dependently by iTachinami et al.ll2011f) . On planets with 
a larger mass than this limit the dynamo shuts down be- 
fore the inner core starts to grow, either because iron 



snow is formed at the core mantle boundary (CMB) or 
the temperature contrast through the CMB falls below 
the value requir ed to ensure convec tion. 

The model bv lGaidos et al.l[2010l is restricted to a dis- 
crete set of masses between 1 and 4.8 M^. In order 
to interpolate the thermal evolution properties to other 
planetary masses we have developed analytical fits for 
Ric, Qconv and idyn- Wc havc found that i?ic and Qconv 
could be fitted using a combination of double exponen- 
tial functions f{x) ^ exp(a;"). It has been observed that 
other complex thermodynamic phenomena could be fit- 
ted well by double exponential functions. This could 
imply that our analytical fits more than having a mere 
practical value could also have some phenomenological 
roots. 

The rising of R-^c for planets below the critical mass 
Merit = 2 A'Iq could be computed analytically using: 



X(i) =Xco (l-exp[-(i-iie)i/'*]) 



(7) 



Here Xoo is a fitting parameter that depends on plan- 
etary mass. The exponents 1/4 and 3 are found almost 
independently of planetary mass. We have additionally 
found that Xoo scales with planetary mass Mp follow a 
simple power law, 



Xoo (Mp) = 1.008 (Mp/Me) 



-7/5 



(8) 



where we have approximated the exponent to the near- 
est rational value (-7/5). Using eq. ^ and ([8]) we can 
calculate analytical approximations of the inner core ra- 
dius for any planetary mass and at any given time after 
the inner core nucleation (solid lines in the upper panel 
of figure [1]) . As expected the combined fit of Xoo and 
x{t) have some discrepancies with the numerical results. 
However, we have verified that most of them come from 
the difference between the actual value of the time of 
inner-core nucletion ije and the fitted value of this quan- 
tity (red line in the lower panel of figure [1} . The observed 
differences, as measured on the time axis, are no larger 
than 0.1 Gyr, i.e. they are below the resolution of the 
magnetosphere evolution model presented here. 

Qeonv is more complex. We have fitted the decaying 
initial phase (release of sensible heat) and the rebounce 
phase after inner core nucleation (release of combined 
sensible and latent heat) independently. The resulting 
piecewise fitting function reads as: 



1012 



9S1 exp(t/tsi)\ 
9L (t-iie)'/^x 

X exp[(i - iic)/iL 
9Sh exp[(t - toh)/ish] 



1/4 



M < Merit, t < U. 

M < Merit, i > U. 
M > Merit, 



(9) 

The fittmg parameters gsi, qsh, Ql, ^si, ^l, ioh and ^sh 
scale also with planetary mass and the scaling coefficients 
are presented in table [1] 

The combination of analytical fits and power-law scal- 
ing of the fitting parameters allows us also to extrapolate 
the thermal evolut ion results to other planetary masses. 
iGaidos et al.ll2010l used their model to study the case of 
planets with a larger mass than 4.8 focusing on the 
dyna mo lifetime. On the other hand. ITachinami et al.l 
12011! also applied their model for planets with a smaller 
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X 






Mass Range 








Me 




Planetary 


Propert: 


les 


Hp 


6378 km 


0.265 


0.7-10 


Rc 


1370 km 


0.243 




Pc 


11 g/cm"^ 


0.271 




Hic 


tic 


2.620 Gyr 


1.360 




Xoc 


1.008 


-0.140 




Q?conv 


gsi 


1.090 


1.528 


0.7-Mcrit 


tsi 


2.907 Gyr 


-0.155 




Qh 


13.529 


1.170 


55 


tL 


0.190 Gyr 


-0.883 




9Sh 


12.616 


-0.350 


Mcrit-10 




-46.096 Gyr 


-2.272 






15.105 Gyr 


-0.965 




TABLE 1 



Scaling coefficients for planetary properties and thermal 
evolution fitting parameters. property x is scaled using a 

SIMPLE POWER LAW FUNCTION X = a^^Mp / Mf^)!^" WHERE Clx IS 
THE REFERENCE VALUE FOR X (THE VALUE AT Mp = 1 Me) AND 
IS THE SCALING EXPONENT. M„it IS THE MAXIMUM MASS FOR 
PLANETS TO DEVELOP A SOLID INNER CORE BEFORE THE SHUTTING 
DOWN OF THE CORE DYNAMO. 

mass than the Earth. We assume here that there are no 
other critical phenomena that avoid the extrapolation 
of the behavior observed in the reported mass-range to 
higher and lower planetary masses. We will model here 
planets in the mass range 0.7 — lOM^. 

Other application of the analytical fitting functions to 
the thermal evolution of TPs is that it provides a way 
to test the impact that different parameters of the ther- 
mal evolution have on the protective properties of the 
evolving PMF. We can check, for example, what would 
happen if an improved thermal evolution model predicted 
a different critical mass Merit or different values for the 
time of inner core nucleation tjc. We will return to this 
sensitivity check in section [5] 

2.2. Planetary magnetic field 

In recent years improved numerical experiments have 
constrained the full set of possible scaling laws used to 
predict the properties of planetary and ste llar convection- 
driven dynamos (see IChristensenI 120101 and references 
therein) . It has been found that in a wide range of phys- 
ical conditions the global properties of convection-driven 
dynamos can be expressed in terms of simple power-law 
functions of the total convective power Qconv available 
for dynamo action. 

One of the most important results of power-based 
scaling laws is the fact that the volume averaged mag- 
netic field intensity inside the convecting shell i?rms = 
(1/ y) f B'^dV does not depen d on the rotation rate (eq. 
6 in lZuluaga fc Cuartasl[20T2l) . 

Brms « CBnns tl^^Ti^} ' \D /Vf QIH, (10) 

where Cerms — 0.24 for dipolar dominated dynamos 
and Cerms = 0.18 for multipolar dynamos, "p^, D ^ 
Rc{l — x) and V « 4/37ri?^(l — x^) are the average 
density, height and volume of the convective shell. We 
are assuming that the whole external liquid iron core is 
convecting. In a real case only a fraction of the core 
volume is involved in dynamo action and therefore the 
magnetic fields predicted with equation [TOl and with our 



assumption will u nderestimate the actual field strenght 
(|Gaidos. E.I 120111 ). We have, however, verified that this 
effect is only important for time periods much longer 
than the time it takes to start the inner core nucleation. 
For a planet with the same mass as the Earth the time 
during which the stellar aggression is the largest is much 
less than that time for inner core nucleation (see section 

a. 

The dipolar field intensity at the planetary surface, 
and hence the dipole moment of the PMF, can be esti- 
mated if we have information about the power spectrum 
of the magnetic field at the core surface. Although we 
cannot predict the relative contribution of each mode 
to the total core field strength, numerical dynamos ex- 
hibit an interesting property: there is a scalable adimen- 
sional quantity, the local Rossby number Ro^ , that could 
be used to distinguish dipolar dominated from multipo- 
lar dynamos. The scalin g relation for i?Oj* is (eq. 5 in 
iZuluaga fc Cuartasll20ll ): 

Rot - Cro, p^-'/'R-'/'D-'/'V-'/^Ql^l^py^. (11) 

Here CroI = 0.67 is a fitting constant and P is the 
period of rotation. It has been found that dipolar dom- 
inated fields arise systematically when dynamos have 
Roi < 0.1. Multipolar fields arise in dynamos with val- 
ues of the local Rossby number close to and larger than 
this critical value. From eq. PT|) we see that in general 
fast rotating dynamos (low P) will have dipolar domi- 
nated core fields while slowly rotating ones (large P) will 
produce multipolar fields and hence fields with a much 
lower dipole moment. 

It is important to stress that the almost independence 
of firms on rotation rate, together with the role that ro- 
tation has in the determination of the core field regime, 
implies that even very slowly rotating planets, for exam- 
ple those whose rotation is locked by the action of the 
tidal effect of its host star (tidally locked planets), could 
have a comparable magnetic energy budget to rapidly 
rotating planets with similar size and thermal histories. 
In the former case the magnetic energy will be redis- 
tributed among other multipolar modes rendering the 
core field more complex in space and probably also in 
time. Together all these facts introduce a non-trivial de- 
pendence of dipole moment on rotation rate very differ- 
ent than that obtained with the traditional scaling laws 
use d by previous works (see e .g. iGriefimeier et al.ll2004l 
and IKhodachenko et al.|[2007D . 

Using the value of Brms and i?o^ we can compute the 
maximum dipolar component of the field at core surface. 
For this purpose wc use the maximum dipolarity fraction 
/dip (the ratio of the dipolar component to the total field 
strength at core surface) that for dipolar dominated dy- 
nam os /^"^ » 1.0 and f or multipolar ones /^"p^ ~ 0.35 
(see IZuluaga &: Cuartai 120121 for details) . To connect 
this ratio to the volumetric averaged magnetic field i?rms 
we use the volumetric dipolarity fraction b^ip that it is 
found, as shown by numerical experiments, conveniently 
re lated with the maxium va lue of /dip through eq. (12) 
in lZuluaga fc Cuartasl[20T^ 

ijmin j-maa: — 11/10 /"I 0^ 

f^dip — Cbdip/dip l-L^J 

where Cbdip ~ 2.5 is again a fitting constant. Finally by 
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combining eqs. (|10m2p we can compute an upper bound 
to the dipolar component of the field at the CMB: 

i?*^ ■f--, ~, ■ — Srms — ~ Brms (13) 

c ~ Urmn ^^^^^ ^. ,. ^^^^^ ^ ^ 

"dip ''bdip 

The surface dipolar field strength is estimated using, 

B';''iR,) = Bfp(^^J (14) 

and finally the total dipole moment is calculated using 
eq. (H]) for r = Rp. 

It should be emphasized that the surface magnetic field 
intensity determined using eq. ()14p overestimates the ac- 
tual PMF dipole component. The actual field could be 
much more complex spatially. Using our model we can 
only predict the maximum level of protection a given 
planet could have from a dynamo-generated intrinsic 
PMF. 

We show in figure [5] the result of applying the previ- 
ously described method to calculate the dipole moment 
for TPs. We compare them with the static value of the 
same quantity as comp uted by the rotation dependent 
dynamo scaling law by I Sand 119931 that has been com- 
monly used in previous works. The differences between 
both approaches are significant. Not only the estimated 
values of the dipole moment are quite different but the 
general dependence on planetary mass and rotation pe- 
riod is much more rich and complex. Those differences 
have important and previously unknown consequences in 
the magnetic protection of tidally locked and unlocked 
habitable planets. We will return to this point in section 

H 

2.3. Stellar properties 

Once we have determined the PMF intensity, the other 
element required to estimate the magnetosphere proper- 
ties is the dynamical pressure of the stellar wind. How- 
ever, magnetic protection does not only depend on the 
size of the magnetosphere. The scale height of the atmo- 
sphere should be also estimated and compared with the 
magnetopause distance. Highly irradiated atmospheres, 
such as those of close-in planets at early phases of the 
stellar evolution, could be expanded enough to be ex- 
posed to the direct action of stellar wind. In this case we 
will therfore also need to estimate the level of high energy 
flux at the top of the atmopshere of our habitable plan- 
ets. In this section we describe the models used in this 
work to calculate all the relevant properties concerning 
the interaction between the star and the planetary mag- 
netosphere and atmosphere. 

The prope rties of low-mass stars are still very uncer- 
tain (see e.g. lEngle fc Guinanir2011[ ). However, since our 
magnetic field model is able only to determine the max- 
imum intensity of the PMF we will only be interested in 
limits of the stellar properties providing the lower level 
of "stellar aggression" (the minimum stellar wind pres- 
sure and XUV irradiation) . Combining upper bounds for 
the magnetic properties of the planet and lower bounds 
for the stellar aggression will give us an overestimation 
of the overall magnetic protection of the planet. In this 
way if, under our conservative model, a planet results 
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Fig. 2. — Maximum dipole moment predicted with the power- 
based scaUng laws used in this work (u pper and m iddle panels) and 
the rotation dependent scaling law bv 1Sandll993l (lower panel). 

endangered or unprotected, the actual case will be even 
worse. 

2.3.1. Habitable zones and tidally locking limits 

The effects a star has on the planetary environment 
depend on three basic factors: 1) the fundamental prop- 
erties of the star, luminosity L^,, effective temperature 
and radius i?*, 2) the average distance of the planet 
from the star (distance to the HZ) and 3) the relation 
between the rotation and orbital period, in particular if 
they are conmensurable (tidally locked planet) or notSD 
(unlocked planet). 

For the stellar properties w e are using the results ob- 
tained bv lBaraffe et al.lll998l (thereafter BCAH98) who 
calculated the evolution of low- mass stars in a wide range 
of mctallicities including the solar value. In all cases we 
assume the same metallicity for the stars as the Sun. 
For simplicity we are also assuming that the basic stel- 
lar properties, temperature, luminosity and radius, are 
constant during the range of time in which we are cal- 
culating the magnetosphere properties t = 0.7 — 3 Gyr 
(see section [2.3.2p . To be consistent with the purpose of 
computing the best case for magnetic protection, we take 
the stellar properties predicted by the model at the end 
of the time interval, i.e. t = 3.0 Gyr. At that time the 
HZ will be placed at the greatest distance from the star 
especially for stars in the mass range AI^, = 0.8 — 1.0 and 
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therefore the stellar wind pressure and XUV irradiation 
will be minimum. 

To determine the dis tance to the HZ w e use here 
the results obtained by iKasting et al.l 119931 Although 
the detailed atmospheric models bv iKasting 
were calculated for only three s tellar masses, we use the 
parabolic fitting developed by iSelsis et all 120071 to cal- 
culate the inner Oin and outer edges Oout of the HZ of 
arbitrary GKM stars: 



flout 



ainAr-A„Ar2)/L7I^ 
- aout AT - Pout AT^)VlJl'c: 



(15) 



where ai„Q and aoutQ are the inner and outer edges 
of the HZ for the present sun, AT = T — Tq and a 
and /3 are fitting constants. To define our HZ we use 
the conservative limits giv en by the criteria o f "recent 
Venus" and "early Mars" (jKasting et al.lll993[ ). In this 
case a„Q = 0.72, ain = 2.7619 x 10-^ An = 3.8095 x 
10-9, a,:„,fT, = 1.77. a^„t = 1.3786 x lO'^, I3,ut = 1-4286 x 
10-9 (ISelsis et al.ll200l . 

Planets at close-in orbits are affected by the tidal in- 
teraction with the host star. This interaction dampens 
the primordial rotation of the planets leaving them in a 
resonant rotational state where the period of rotation P 
becomes commensurable with the orbital period P^, 

P : Po^n:2 (16) 

Where n is an integer larger than or equal to 2. The 
value of n is determined by multiple dynamical factors 
one the most important being the orbital eccentricity. 
The maximum distance atid at which a solid planet in a 
circular orbit becomes tidally locked before a given time 
t is given by (jPealail977i) . 




D.O 0.4 0.8 

LfL,.,R/R,.,T/T 



a(AU) 



Fig. 3. — Stellar properties, HZ and tidally locking limits for 
solar metallicity GKM stars. Stellar properties are computed from 
the BCAH98 model at t = 3 Gyr. Tidally locking limits assume 
Ppiim = 17 hours, and Q = 100. 

Given the dependence of the magnetosphere properties 
on the stellar wind pressure (equations [3][6]) the rotation- 
independent model will give us values for the magne- 
topause fields, standoff distances and polar cap areas, 
between 10-40% off the values obtained with a more de- 
tailed model (e.g. the extended model bv I Weber fc DavisI 
Il967f ). Magnetopause fields have the largest uncertain- 
ties 40%) while standoff distances and polar cap ar- 
eas, that as we will also shown are the most critical pa- 
rameters, are off by just ~ 10% (standoff radius will be 
overestimated while polar cap areas are underestimated). 

According to the Parker's model the stellar wind aver- 
age particle velocity v at distance d from the host star is 
obtained by solving the Parker's wmd equation: 



atid(0 =0.5AU 



(M,/M0)^Ppri.n ' 

Q 



1/6 



1-1/6 



(17) 



Where the primordial period of rotation Pprim should 
be expressed in hours, t in Gyr and Q is an adimensional 
dissipation function. For the purposes of this work we 
assume the same value of the estimated primordial rota- 
tion period of the Earth for all planet s, Pprim ~ 17 hours 
(jVarga et al.lll998t IDenis et al.ll2011[ ) and a dissipation 
function Q = 100 typical for terrestrial planets. 

In figure[3]we depict the value of the stellar properties, 
HZ and tidally locked limits for the stars in the mass- 
range studied in this work. 

2.3.2. Stellar wind 

The properties of the stellar wind change in time and 
vary with the distance from the host star. Dur to this 
dependence on distance and for the sake of simplicity 
we have used here t he pure hydrodynamical isothermal 
model developed bv lParkeri ri958 (hereafter the Parker's 
model). It has been shown that this simple model re- 
liably predicts the stellar wind properties of stars with 
periods of rotation of the sa me order as the prese nt so- 
lar value, i.e. P ~ 30 days (jPreusse et al.l [20051) . For 
rapidly rotating stars, i.e. young stars and/or active 
dM stars, the Parker's model underestimates the stellar 
wind properties by a factor up to 2 ((Preusse et al..,2005l ). 



log u = 4 log p 



(18) 



where u = v/vc and p 
velocity and distance and v, 



d/dc are a normalized 



yJkBT/m and dc — 
GMi,m/{4kBT) are respectively the local sound veloc- 
ity and critical distance where the stellar wind becomes 
subsonic. In this model T is the temperature of the stel- 
lar corona and appears here as the key parameter that 
determines the properties of the stellar wind. The num- 
ber density n is determined from the velocity at each 
distance using the continuity equation. 



vm 



(19) 



The time dependence of the stellar wind properties is 
much harder t o estimate. We have us ed here the formu- 
las derived bv iGriefimeier et al.ll2004l and lLammer et al.l 
[2004 and that were originally based on the o bservational 
estimates of the stellar mass-loss rate by iWood et al. 
I 2002i and the theoretical evolution models bv INewkirk 
119801 For main-sequence stars and times t > 0.7 Gyr, 
the long-term averaged stellar wind velocity and number 
density as measured at a distance of 1 AU is estimated 
by: 



VlA\j{t) = Wo ( 1 + ^ 



(20) 
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niAu(i) =no[l + 









- - tV.M ' ' D.V;U, ■ ■ 1.0.1/. 











Where a„ = -0.43, 



a„ = —1.86 ± 0.6 and r — 

25.6 Myr ()Griefinieier et al.ll2009l) . vo = 3971 km/s and 
riQ — 1.04 X 10^"'^ m~'^ were estimated assuming present 
long-term averages of the solar wind as measured at 
the distance of the Earth, i.e. n(4.6 Gyr, 1 M©, 1 AU) = 
6.59 X lO^m-s and t;(4.6 Gyr, 1 Mq, 1 AU) = 425 km/s 
(|Schwenn|[l990[ ). For times t < 0.7 Gyr stellar wind 
models are too uncertain and the est imates provided b y 
eqs. ^ and (HU become unreliable (|Wood et al.ll2002[ ). 

Gricfimcicr et al. 2007 devised a clever way to combine 
the distance dependent estimation of the stellar wind 
properties, given for example by the Parker's model, with 
the time variation of the reference number density and 
velocity given by eqs. (pOt and (PTt . For the sake of 
co mpleteness in this work we summarize the procedure 
bv lGriefi meier et aLll2007l but for further details we refer 
to section 2.4 of that work. 

The stellar wind properties at time t and distance d 
for a given stellar mass are calculated by estimating first 
the coronal temperature T{t) for which the velocity ob- 
tained with the Parker's equation evaluated at d = 1 
AU coincides with the reference velocity obtained by eq. 
(pn|) evaluated at time t. Using the coronal temperature 
T{t) the Parker's equation provides the velocity at any 
distance from the star at that time. The number density 
is obtained from the continuity equation (|19p assuming 
that the stellar mass-loss rate scales with the stellar ra- 
dius as Mi,{t) = MQ{t){Ri, / RqY , where the solar refer- 
ence mass-loss rate Mq is computed using: 

MQ(i) =47r(lAU)2 mniAu(t)viAu(t) (22) 

The values of the stellar wind dynamical pressure 
-Pdyn = mnv^ inside the HZ of four different stars as 
computed using the procedure described before are plot- 
ted in the upper panel of figure |4l 

2.3.3. Stellar XUV fluxes 

The XUV luminosity of a star depends on its level of 
chromospheric and coronal activity which in turn depend 
on the rotation rate of the star. It is now well known that 
the rotation of main sequence stars slows down with age. 
It follows from there that the XUV luminosity should also 
decrease monotonically in time. Estimating the variation 
in time of the XUV luminosity for stars of different spec- 
tral types is harder than thought. There is a large disper- 
sion of rotation periods and hence in the XUV lu minosity 
of stars of the same spectral type (jPizzolato et al.ll2003l ). 
Differences in one order of magnitude have been observe d 
in stars with the same age (see e.g. iMicela et al.l 119961 ) . 
Additionally EUV radiation is absorbed by the interstel- 
lar medium and therefore we depend on proxies to esti- 
mate the total XUV luminosity. X-ray luminosity being 
the most common. 

Using observational data and independent estimations 
of stellar ages several authors have developed empiri- 
cal laws providing the value as a function of time of 
different XUV luminosity proxies, e.g. Lx (see e.g . 
iRibas et al.|[2005l iPenz fc Miccla.2008, .Penz et al.ll2"008l 
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Fig. 4. — Evolution of the stellar wind dynamical pressure (upper 
panel) and the XUV flux (lower panel) inside the HZ of a selected 
set of stellar masses. 



iLammer et al.l 120091 iGarces et al.l 1201 ID . Given the im- 
plicit uncertainties in the estimation of ixuv and to be 
consistent with our goal of obtaining the best case for 
magnetic protection, we have selected the model pre- 
dicting the lowest values of the XUV lumino sities. The 
empirical law obtained bv IGarces et al.|[2011l is the best 
suited for that purpose. According to that result the X- 
ray luminosity of GKM stars change over time following 
the simple power-law function: 



L 



X 



6.3 X lO-^L* if t < ti 
1.89 X 10^8 otherwise 



(23) 



where L^, is the bolometric luminosity of the star and 
ti is the end of the saturation phase that scales with i^, 
as. 



ti = 0.06 Gyr 



Lq 



-0.65 



(24) 



Luminosities are in units of erg s ^ . For simplicity we 
assume £xuv ~ -^x- Using this model the Present Earth 
Value of the XUV flux (thereafter PEL) is 0.64 erg cm'^ 
s^^ which, as expecte d, underestimates the observed 
value of this quantity (jJudge et al.l l2003t iGuinan et al.l 
12001 . 



Evolution of magnetic protection in habitable terrestrial planets 



9 



In figure 13] we plot the value of the XUV flux as a func- 
tion of time as measured in the HZ of four different stars. 
XUV fluxes in the range of 10 - 1000 PEL seem to be 
common in early phases of stellar evolution posing severe 
constraints on the survival of planetary atmospheres of 
unmagnetized and weakly magnetized planets. We will 
return to this point in sectional 

3. ATMOSPHERIC THERMAL EXPANSION AND 
MASS-LOSS 

One of the most critical points regarding the survival 
of the atmospheres of low-mass planets is the fact that in 
the early phases of stellar evolution they would be under 
extreme conditions of X and EUV irradiation. A highly 
irradiated and relatively light atmosphere will have high 
levels of thermal and non-thermal mass-losses especially 
if it is unprotected against the action of the stellar wind. 
It has been estimated that under no or even a weak pro- 
tection from an intrinsic magnetic field a TP could lose 
its atmosphere or most of its volatile content in a time- 
scale much sh orter than that requir ed for the evolution of 
life (see e.g. (jLammer et al.|[2012D ). Understanding the 
relation between the XUV-induced atmospheric expan- 
sion and mass-loss, and the properties of the planetary 
magnetosphere is of fundamental importance in assessing 
the problem of magnetic protection. 

The key property in distinguishing between a mag- 
netically protected atmosphere and an exposed one is 
the radius of the exobase i?oxo, defined as the distance 
where the mean-free path of atmospheric particles could 
be comparable to the size of the planet. This is the 
limit where atmospheric particles, given the proper en- 
ergetic or flux conditions, could escape from the plan- 
etary atmosphere. The height of the exobase depends 
on many complex factors ranging from the opacity of 
the atmopsheric gasses to the high energy radiation from 
the star, the interaction between the charged and neu- 
tral components of the high atmosphere and a complex 
network of chemical and photochemical reactions of the 
atmospheric con stituents (for a complete description see 
iTian et al1l2008( ). All these factors are critically deter- 
mined by the chemical composition of the atmosphere. 
In the last few years several authors have, using de- 
tailed chemical, thermal and hydrodynamical models of 
TPs atmospheres, calculated the exospherc properties for 
two chemical compositions: N2 rich or Earth-like com- 
posit i on atmospheres fiWatson et al. 1981; Kulik ov et al.l 
l2006t ITian et al.l [ 200l and drv CO2 or Venus-like com - 
position atmospheres (|Tianl [20091 : iLammer et al.ll2012f ). 

In order to include this important component in our 
comprehensive model we have used here the exospheric 
properties computed with two diffe rent models: a N? rich 
atmosphere of an Earth-like planet (|Tian "etaII[2008P and 
a CO 2 rich atmosphere of massive SEs Mp > (jTianI 
|2009() . The radius of the exobase in these models as a 
function of the XUV flux is depicted in figure [S] We see 
that N2 rich atmospheres of Earth-like planets expand 
further than ~ 10i?p at XUV fluxes larger than ~ 20 
PEL. In this case the atmosphere become exposed even 
if they have the maximum level of magnetic protection, 
i.e. Rs ^ IQRp and could be stripped off by thermal 
and (stellar-wind related) non-thermal mass-losses. On 
the other hand, CO2 rich atmospheres are cooled off ef- 
ficiently through the emission of IR-radiation in the 15 




Fig. 5. — Radius of the exobase for the models used in this work. 
Left panel corresponds to an Earth-mass planet with a N/O-rich 
atmosphere. Right panel shows the exobase radius for super-Earths 
with C02-rich atmospheres. Dotted lines in the right panel are the 
linearly interpolated and extrapolated exobase radius for masses 
from top to bottom 3.5, 4.5, 5.5, 7.0 and 9.0 M®. XUV fluxes are 
expressed in terms of the present Earth level, PEL 5 mW/m^. 

/im CO2 band and become able to withstand higher XUV 
levels. 

In order to evaluate if a TP is magnetically protected 
we will compare the evolving value of the exobase radius 
and the standoff distance computed with our dynamical 
magnetosphere model. It is expected that at the earli- 
est phases of stellar and planetary evolution the exobase 
radius would be larger than the standoff distance irre- 
spective of the existence of an intense early dynamo. To 
expect that the magnetosphere protects the planetary at- 
mosphere from the very beginning is simply unrealistic. 
Therefore the critical property to evaluate the level of 
magnetic protection of a given planet is the time A^oxp 
during which this exposition state is maintained. More 
precisely the total atmospheric mass lost during this in- 
terval Moxp will give us the information we require to 
determine if the planetary atmosphere can survive the 
early "aggression" of its host star. We return to these 
important properties of the magnetic protection evolu- 
tion in the next section. 

4. RESULTS 

Combining all the elements of our comprenhensive 
model we have calculated the evolving magnetic protec- 
tion conditions of potentially habitable TPs. 

We have caclculated and compared the level of mag- 
netic protection for planets in the mass range 0.7— 10 
in three different cases: 1) instantaneous values of the 
magnetosphere and exosphere properties, especially at 
the earliest phases of the thermal evolution (t w 1 Gyr), 
2) evolution in time of the same properties during a 
time-scale comparable to the development of complex 
life (i = 1 — 3 Gyr) and 3) the cummulative atmospheric 
mass-loss of major constituents in the same period of 
case 2. 

We study these three cases for three type of planets: 

a) hypothetical planets in the mass-range 0.7 — 10 M^, 

b) an Earth-twin, i.e. a habitable planet with the same 
mass and period of rotation than the Earth, but orbiting 
different type of stars and c) the Earth-like planets and 
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SEs already discovered inside the HZ of their respective 
host stars, including, for reference purposes, Venus and 
the Earth itself. In table [2] we summarize the properties 
of the planets in the last group. 

To include the effect of rotation in the properties of 
the PMF we have assumed that planets in the HZ of 
late K and dM stars (M < O.QMq) are tidally locked 
and therefore their periods of rotation are equal to their 
orbital periods {n—2 in eq. [T6| . For planets that preserve 
their primordial periods of rotation wc assume values for 
P in the range 1 — 100 days as predicted by models of 
planetary formation (Migu el fc Bruninill2010( ) with P = 
1 day as the preferred reference value. 

Figures [6] and [7] show the evolution of magnetosphere 
properties for tidally locked and unlocked habitable plan- 
ets. In all cases we have assumed that the planets are in 
the middel of the HZ. Even in early times tidally locked 
planets of arbitrary masses have a non-negligible magne- 
tosphere {Rs > 1.5 Rp). As the star evolves the dynam- 
ical pressure of the stellar wind decreases more rapidly 
than the dipole moment (see figures [2] and |4]) and the 
standoff distance grows. The critical boundary observed 
in the middle and rightmost panels are a product of the 
inner core nucleation. Planets to the right of the bound- 
ary created by a concentration of isolines still have a 
liquid core and therefore develop lower dipole moments. 
On the other hand, the inner core in planets to the left of 
the same boundary have already started to grow before 
that time and therefore their dipole moments are much 
larger. 

Previous estimates of the standoff radii for tidally 
locked plane ts are lower than the valu es reported here. 
For example. iKhodachenko et al]|2007l place the range of 
standoff distances well below 2 Rp even under mild stel- 
lar wind conditions (see figure 4 in their work). This is 
easily explained since they also underestimate the maxi- 
mum dipole moment for this type of planet. While they 
predict maximum dipole moments for tidally locked plan- 
ets around stars with M^, « 0.5Mq in the range of 0.022 
to 0.15 A^e, our model predicts dipole moments as large 
as 0.8 A^e for planets with the largest mass and at times 
as early as t ss 1 Gyr. 

A lower standoff distance means a larger polar cap. 
Tidally locked planets, even under the assumption that 
their exobases are not larger than the magnetosphere, 
have more than 15% of their atmospheric surface area ex- 
posed to open field lines where thermal and non-thermal 
processes could efficiently remove atmospheric gasses. 
Moreover, our model predicts that this type of planet will 
probably have multipolar PMF (thereafter paleomagne- 
tospheres) which only increase the areas where field lines 
are op en to the interplaiietary space and magnetotail re- 
gions fStadel mann et aT]l2010t ). 

Unlocked planets (figureEj) as expected seem to be best 
protected by extended magnetospheres Rs > ^ Rp and 
lower polar cap areas < 10%. It is interesting to 
notice that in both cases and at early times t ~ 1 Gyr 
a smaller planetary mass also implies a lower level of 
protection. This fact seems to contradict the idea that 
low-mass pla nets {Mp < 2.5) are best suited to develop 
intense PMF liGaidos et al.ll2010HTachinami et al.ll2011l: 
iZuluaga fc Cuartasll2012D . This apparent contradiction 
is explained by taking into account that the early dy- 
namo generated magnetic fields of planets with very dif- 



fer en^jnassesareof the same order (see e.g. figure 8 
in lGaidos et al.|[2010D : although larger planets produce 
more convective energy, the oversized liquid core and a 
larger planetary radius produce a surface magnetic field 
similar in intensity to that of a smaller planet. Hav- 
ing similar surface PMF, planets with larger masses will 
have much larger dipole moments, M ^ Rp and under 
the same stellar conditions will be best protected against 
the action of the stellar wind. This situation is not longer 
valid when low-mass SEs (Mp < 2.5) develop an inner 
core and the available energy for convection is largely 
increased. 

We have also calculated the standoff distance and polar 
cap area for different evolutionary stages of the already 
discovered habitable planets enumerated in table [2] The 
position are indicated with circles, whose size is propor- 
tional to their measured or estimated planetary radii, 
in the contour plots of figures [6] and [7l Planets whose 
properties are out of the ranges used in these figures are 
indicated with arrows at the border of each panel. 

A more detailed account of the evolution of magneto- 
sphere properties for already discovered habitable plan- 
ets is presented in figure [H The values of magnetosphere 
properties at the present age of the planet are indicated 
with circles whose size is proportional to the planetary 
radii. The case of Venus and that of the Earth are also 
included for reference purposes. 

In all cases we have used values for the period of ro- 
tation ranging from a minimum, corresponding to an as- 
sumed primordial rotation rate of Pmin = 24 h (with 
exception of the Earth for which we have used Pmin = 17 
h), and a maximum, equal to the orbital period (tidally 
locking) . The minimum and maximum values of the ro- 
tation period determine the upper (lower) and lower (up- 
per) boundaries of the shaded regions in the upper panel 
(lower panel) respectively. 

The case of Venus is particularly interesting in order 
to analyse the rest of the planets. The dynamo of Venus 
probably shut down at t = 3 Gyr as a c onsequence of the 
drying of the mantle (jChristensen et a l. 2009). A mas- 
sive loss of water induced by a runaway greenhouse and 
unsufficient early magnetic protection played a central 
role in the extinction of the early Venusian PMF. The 
case of Gl 581 d, GJ 667Cc and HD 85512b, though sim- 
ilar to Venus, is much more complex. On one hand their 
masses are larger than Venus' and therefore their gravity 
could provide additional protection to the atmospheric 
mass-loss. On the other hand, their dynamos shut down 
in times t ss 3 — 5 Gyr exposing them to the direct ac- 
tion of the stellar wind. In G and K stars this situation 
could not be consider a big threat since the stellar wind 
had also decreased its intensity when the planet lost its 
PMF (see figure [4]), this will also be the case for HD 
85512b. However, Gl 581d and GJ 667Cc are located 
at the HZ of dM stars where the stellar wind pressures, 
even at times as late as 4 Gyr, are intense enough to 
massively erode their atmospheres or to make them lost 
their volatile content. 

As explained in section [3] the magnetic protection of 
a habitable planet is not only a function of the magne- 
topshere properties. In order to asses the problem we 
need also to evaluate the effect of the XUV radiation in 
the outer atmosphere expansion. In figure [5] we compare 
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Planet 






a(AU) 


Fo (days) 


e 


S-type 






age(Gyr) 


tid. locked 


Earth 


1.0 


1.0 


1.0 


365.25 


0.016 


G2V 


1.0 


1.0 


4.56 


No 


Venus 


0.814 


0.949 


0.723 


224.7 


0.007 


G2V 


1.0 


1.0 


4.56 


Probably 


GJ 667Cc 


4.545 


1.5* 


0.123 


28.155 


< 0.27 


Ml. 25V 


0.37 


0.42 


> 2.0 


Yes 


HD 85512b 


3.496 


1.4* 


0.26 


58.43 


0.11 


K5V 


0.69 


0.53 


5.6 - 8.0 


Yes 


Gl 581d 


6.038 


1.6* 


0.22 


66.64 


0.25 


M3V 


0.31 


0.29 


4.3 - 8.0 


Yes 



TABLE 2 

Properties of already discovered SEs inside the HZ of their host stars. For reference purposes the properties of Venus 
AND the Earth are also included. Values of radii marked with an * are unknown and were es timated using the mas s-radius 

RELATION FOR PLANETS WITH THE SAME COMPOSITION AS THE EARTH, I.E. Rp = R^{Mp / M^)'^-'^'' HVALENCIA ET AL.II2006I) . 



1,0 Gyre f= 3.0 Gyre / = 4.5 Gyre 




(-1/.) (.1/.) (.U.) 



Fig. 6. — Evolution of the standoff radius and polar cap area for tidally locked planets: the period of rotation of each planet is assumed 
equal to the orbital period at the middle of the HZ of low mass stars Mi, < 0.6Mq (ticks on the right vertical axis). The position of already 
discovered potentially habitable TPs is indicated with a circle or with an arrow when their properties fall beyond the axis limits. 



the evolution of the exobase radius and the standoff dis- 
tance for two representative TPs planets: an Earth-mass 
planet with a N/O-rich atmosphere in the HZ of G-K 
stars and a super-Earth with Mp = and a C02-rich 
atmosphere in the HZ of dM stars Mi, < 0.6. For the 
Earth-mass planet we have assumed periods of rotation 
in the range P — 1 — 1 days which are compat ible with 
TP formation theories (jMiguel fc Bruninill2010l ). For the 
rotation rate of the super-Earth we have assumed rota- 
tion periods of between 1 day (primordial rotation) and 
the orbital period (tidally locked case). Both planets 
were placed at the middle of the HZ of each star. Stand- 
off distances for the HZ of stars in the mass-range studied 
in each case do not change significatively. For reference 
purposes we have only plotted this quantity for the star 
with the largest mass in the interval considered in each 
case (1.0 Mq for the upper panel and 0.5 Af© for the 
lower panel) 

For a solar-mass star and P = 1 day (dipolar dom- 
inant early PMF), the exposure time is Atcxp ~ 200 
Myr. However, if the primordial period of rotation is 3-5 
days (multipolar field) the exposure could be increased 
by up to 1 Gyr threatening the atmospheric stability or 
its content of volatiles. The exposure time for an Earth- 
mass planet increases monotonically when it is located at 
closer distances to late G and early K stars. If the Earth- 
mass planet is located in the middle of the HZ of a star 
with Mi, < 0.86 (oHz ~ 0.80 AU) it would be subject to 



levels of XUV irradiation high enough to blowing off the 
atmosphere. This fact points to the existence of a sub- 
region of the HZ that we can call a Magnetic-restricted 
Habitable Zone (MHZ), where planets with a given mass 
and atmospheric composition could preserve their atmo- 
spheres. It is interesting to note that although Venus 
was inside the solar HZ during the first Gyr of the solar 
system evolution it has always been outside of the MHZ 
for Earth-like planets with N/O-rich atmospheres. CO2 
rich atmospheres which are less propense to expand, i.e. 
i?cxo < 2 for Fxvv ~ 100 (Kulikov et al. 2006) are prob- 
ably protected at the critical early phases of planetary 
evolution, provided the planet preserves the conditions 
required for dynamo action (e.g. an hydrated mantle, 
mobile lids, etc.) In that case the MHZ will extend to 
stars with lower masses or closer habitable distances. 

The case for massive super-Earths in the HZ of dM 
stars has several interesting differences with respect to 
the magnetic protection of Earth-mass planets. The con- 
tinued activity of very late dM stars (M^ ~ 0.1) will 
produce a constant level of XUV irradiation and hence 
a constant exobase radius during times longer that the 
dynamo life-time itself (magenta solid line). Under this 
condition and assuming that the planet is tidally locked 
in the first Gyr, the atmosphere will be permanently ex- 
posed and will probably be subject to large thermal and 
non-thermal mass-losses. However, for stars with masses 
Mi, > 0.2 a tidally locked planet will only be exposed 
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Fig. 7. — Same as figure[6]but for unlocked planets with constant period of rotation P = 1 day. 



Earth ( 1.0 ;!/ ) — GJ 667C c (4.5 ;i; ) — Gl 581d (6.0 ,» I 
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Fig. 8. — Evolution of magnetosphere properties for already dis- 
covered habitable SEs. Shaded regions are limited by the properties 
determined with the minimum (P fa 1 day) and maximum period 
of rotation of the planet (P Po). 



durin g few Myr t o 1.0 Gyr. If the mass-loss is not too 
large (iTianI [20091) the planetary atmosphere could sur- 
vive the initial stellar aggression. Again we could iden- 
tify here a minimum stellar mass (HZ distance) for which 
the atmosphere is critically exposed. For a tidally locked 
6 super-Earth that minimum mass is close to 0.3 Mq 
(anz ~ 0.15 AU). The inner limit of the so-called MHZ 
is farther away for SEs with larger masses but still covers 
a significant fraction of the tidally locked region of the 
HZ. 

In order to quantify in more detail the effect that the 
early exposition to the stellar wind has on the atmo- 
sphere erosion of habitable TPs we have calculated the 
exposure time and the total mass loss during this crit- 
ical period for an Earth-mass planet in the HZ of G-K 
stars and tidally locked SEs. The results are presented 
in figure [TUl 

In the case of super-Earths we have use d the flux of 
carbon reported in figure 8 of [Tia^[20Ql For masses 
not included in the simulations we have linearly interpo- 
lated and extrapolated the already reported results for 
planets with masses of 6, 7.5 and 10 M(q. The flux of 
oxygen, the most abundant ion in higly irradiatend N/0- 
rich atmospheres, in the case of the Earth-mass planet, 
was calculated multiplying the estimated number den- 
sit y n as estimated fr om the exobase definition (eq. 18 
in kulikov et al.ll200C) and the exobas e bulk velocity re- 
ported in figure 8d of lTian et al.l l2008 In both cases the 
total mass-loss is compared with the assumed initia l con- 
tent of CO2 (2.5 X lO''^ for an Earth -mass planet, iTianI 
|2009() and the amount of water in the planetary hydro- 
sphere (4.6 x 10'*^ molecules for an Earth-mass planet, 
IKulikov et aLll2006D . For planets with a mass larger than 
that of the Earth and for the sake of simplicity these ini- 
tial content of critical volatiles were scaled linearly with 
mass. 

We have observed that exposure time and mass-loss, 
depend on the distance from the star but are almost in- 
dependent of the stellar mass compatible with that dis- 
tance. Therefore, an Earth-mass planet placed at ^ 0.8 
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Fig. 9. — Radius of the exobase (solid lines) and standoff dis- 
tance (shaded area) as a function of time for two potentially hab- 
itable TPs. Upper panel: Earth-mass planet with a N/O-rich at- 
mosphere. The shaded region corresponds to the standoff distance 
computed for a solar-mass star and assuming periods of rotation 
of 1 day (upper limit) and 10 days (lower limit). Dashed line is for 
P = 1.5 days and dash-dotted line is for P = 3 days. Lower panel: 
massive super-Earth with a C02-rich atmosphere. The shaded re- 
gion corresponds to the standoff distances calculated in the HZ of a 
star with O.SMq and periods between 1 day (upper limit) and the 
orbital period (tidally locked case, lower limit). Dashed line corre- 
sponds to a rotation period P = 4 days or approximately l/17th 
of the orbital period. 



& 2.0 



10" 



10" 



1.34AU (1.05 M, ) 
I.IBAU (1.00 .U ) 



1.03 AU (0.95 ,U I 
OSIAU (0.90 .U I 



0.81 AU (0.86 




_i_ 



0.04 
0.10 



AU(0.10fl/,) 
AU (0.20 jU .) 



0.15AU(a30,Ul — 0,25 AU (0.47 j;) 
0.20 AU (0.40 M ,) 




1 2 3 4 5 6 7 



of planetary mass for two atmospheric compositions: N/O-rich at- 
mospheres (Earth-mass planet) and C02-rich atmospheres (super- 
Earths). Planets were placed at different distances from their host 
stars (see legends). For each distance the mass of the star, where 
that distance corresponds to the middle of the HZ, is indicated 
in parenthesis. For the Earth-mass planet we assumed a constant 
period of rotation P = 1 day. The super-Earths were placed in the 
HZ of dM stars and are tidally locked, i.e. P = Pq. For reference 
we have included the duration of the dynamo tdyn with respect to 
the initial time (0.7 Gyr) (upper panel, dashed line) and the CO2 
and H2O content scaled with the Earth's mass (dahed and dash- 
dotted lines in the lower panel). The exposure time and mass- loss 
predicted for the already discovered habitable planets are marked 
with squares. 



AU will be subject to a comparable exposure around a 
0.75 Mq where that distance correspond to the outer 
boundary of the HZ (see figure or around a 0.86 Mq 
star (middle of the HZ) or around a 1.05 Mq star (inner 
boundary of the HZ). This is the reason why we have 
decided to parametrize the results in terms of the plane- 
tary distance from the star rather than in terms of stellar 
mass as in figure |9l 

The exposure times for the Earth-mass planet are com- 
patible with the results presented in figure [HI The early 
mass-loss for this type of planet ranges from 25% of the 
scaled mass of the ocean for a distance 20% larger than 
present Earth, to a value larger than this critical thresh- 
old for a planet 20% closer than the Earth. The exposure 
time and mass-loss for tidally locked super-Earths de- 
crease with mass as expected from the decreasing of the 
exobase radius with an increasing planetary mass (see 
figure [5l We confirm here the existence of a minimum 



distance for a given planetary mass beyond which the 
carbon mass-loss is below the expected CO2 content. For 
example, a super-Earth with mass 6Mq located at dis- 
tances larger than 0.1 AU (long dashed lines in the lower 
panel) will always have mass-losses below that threshold. 
This limit corresponds to the middle of the HZ around 
a 0.2AfQ star which confirms the analysis derived from 
figure [5] 

With the information at hand and assuming that the 
already potentially habitable super-Earths have similar 
compositions to the Earth and CO2 rich atmospheres we 
can conclude that GJ 667Cc has already lost a signifi- 
cant fraction of its atmospheric mass and it is probably 
now uninhabitable. Gl 581d and HD 85512b seem to be 
in a safe region of the parameter space. Altough our 
model provides a very conservative estimate of the early 
exposure conditions, the estimated mass-losses for both 
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planets are more than one order of magnitude below the 
critical threshold. Even if we accept that their dynamos 
are weaker or if we include the additional mass-loss pro- 
duced after their dynamos shut down (see figure [8]) the 
total mass-losses are still below the scaled CO2 content. 

5. DISCUSSION AND FURTHER ANALYSIS 

This is the first attempt to integrate into a single com- 
prehensive model all the relevant physical phenomena 
involved in the magnetic protection of habitable planets. 
Although several components of the model are expect- 
ing important improvements in the following years (bulk 
properties of low-mass stars and their evolution, struc- 
ture and evolution of stellar winds, physics of highly ir- 
radiated atmospheres, atmospheric escape in magnetized 
planets) the esential elements have been put together to 
produce an integrated view of the role of magnetic fields 
in the survival of the atmosphere of habitable TPs. 

One important source of uncertainties in our model, 
especially when applied to already discovered TPs, is 
the assumption that all of them have similar composi- 
tions to the Earth. Planets with elemental and miner- 
alogical compositions different to our planet are proba- 
bly in ore abundant than previously thought (jBond et al.l 
l201Clf ). Although numerical models of the bulk interior 
of solid planets with ve ry different compos itions have al- 
ready been computed (jSeager et al.l 120071 ) the detailed 
structure, mincralogical phases, thermal profiles and evo- 
lution, among other geophysical relevant information, are 
waiting to be studied in more detail for this type of 
planet. Plate tectonics, mantle rheology, additional inte- 
rior heating sources (e.g. radioactive and tidaly heating) 
and the formation, composition and thermal structure 
of a metallic core are key properties to study the ther- 
mal and magnetic field evolution on planets with very 
d ifferent composition to the Earth. 

iGaidos et al.l ()2010D have studied the effect that an in- 
creased core size (larger Fe/Si ratio) or a different amount 
of radionuclides have on the thermal evolution and hence 
the PMF evolution of planets with different composition 
than the Earth. They found that increasing the Fe/Si 
ratio for an Earth-mass planet, i.e. increasing the ra- 
dius of the core, will have two effects on the thermal and 
magnetic field evolution (see figure 7 in their paper): 1) 
an earlier formation of the solid iron core and 2) more 
intense surface fields (mainly due to a reduction in the 
distance between the planetary and core surface without 
a significant change in the convective power). A change 
in the predicted magnetic field intensity of almost one 
order of magnitude was observed changing the core ra- 
dius between 0.5 and 1.5 times the Earth's core radius. 
Using equations [T] and [5] it is predicted that a larger rel- 
ative amount of Fe could imply a standoff distance of up 
to 3 times larger than that predicted for a planet with 
a composition similar to Earth. Also an earlier inner 
core formation would imply increased magnetic protec- 
tion during the critical earlier phases of stellar and plan- 
etary evolution. The effect of a different amount of ra- 
dion uclides on the PMF strength is negligible (see figure 
7 in IGaidos et al.ll2010l ). In all these cases the magnetic 
protection conditions predicted here will be improved. 
Thus, for example, GJ 667Cc will have better chances 
of preserving its atmosphere if its content of Fe is much 
larger than that of the Earth. 



The effect of different rheological properties and 
thermal structure on SEs have also been studied in 
the detailed mantle a nd core thermal evolution by 
iTachinami et al.l ()2011| ). Although they also assumed 
the same elemental and mincralogical composition as the 
Earth the attention they put on the role of variations 
in the rheological properties of the mantle and different 
thermal conditions at the core mantle boundary (CMB) 
could be helpful in roughly guessing what could happen 
on planets with different compositions. They have found 
that the total lifetime and surface strength of the PMF 
ar e sensitive to changes i n these properties (see figure 11 
in ITachinami et all 120 111) . Changes become particularly 
important for planets with masses larger than the Earth 
{Mp > 2Mq). For example, they predict that the metal- 
lic core of planets with a mass as large as 5 could 
actually cool enough to develop an inner solid core, pro- 
vided that on one hand the viscosity of the mantle is 
weaker dependent on pressure and on the other hand, 
that the temperature contrast through the CMB is up to 
10 times larger than that expected for the Earth. In this 
case larger-mass-planets could develop strong PMF and 
its magnetic protection could be ensured even at closer 
distances to their host star than those predicted here. 

We have simulated the previously discussed effects by 
changing the fitting parameters of the phenomenological 
functions describing the thermal evolution (equations [7]|9] 
and scaling parameters in table [Ij . We present in figure 
ITT] the results of this "test of sensitivity" to our model 
of variations in the thermal evolution details. We have 
plotted the percentual variation of the standoff distance 
and the factor of variation in the total mass-loss when the 
following properties of the thermal evolution are modi- 
fied with respect the nominal values used in this work: 
1) the radius of the metallic core Rc (scaling parameter 
Q^Rc), 2) The time of inner core nucleation tiC (scaling 
parameter atic) and 3) the critical mass for solid inner 
core nucleation. Merit. 

We see that the value of standoff distance predicted by 
the model is very robust. Most of the modification of the 
model parameters produce changes in the nominal value 
of only ~ 20% or less. This could be explained by the 
weak dependence of the standoff distance on the dipolar 
moment, eq. ([5]) Rs ~ M^^^. Only when the core radius 
is increased (we have used Rc ~ 1.2i?p where Rq is the 
nominal core radius scaled with the mass of the Earth) 
is a difference in the standoff radius of 40% observed in 
the case of massive SEs. For planets with lower masses 
the difference is still around the 20% limit. 

The exposure time and mass-loss are much more sen- 
sitive to changes in the thermal evolution parameters. 
This is explained taking into account that even small 
variations in the standoff distance could increase by a 
factor of 2-5 the time during which the planet is exposed 
directly to the stellar wind erosion. As expected the 
largest changes are observed when we change the core 
radius (Fe/Si ratio). The predicted mass-losses for EPs 
and massive SEs are almost 10 times lower than those 
obtained with the nominal model. SEs with intermedi- 
ate mass {Mp ^ 4M0) are less sensitive to core radius 
changes and the mass losses decrease by a factor of 2 
only, when the core radius is increased. These results 
support to the conclusion that Gl 581d and HD 85512b 
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Fig. 11. — Test of the model's sensitivity to changes in the param- 
eters of the thermal evolution model {Rc, tic a-nd ^crit) ^^nd the 
stellar wind and XUV irradiation properties (coronal mass ejection 
conditions or CME and Fxuv)- The value of the magnetic pro- 
tection key quantities Rg, texp and A^exp are compared with their 



nominal values R^, 



and Af£,„ as calculated with the set of 



exp ""^^^ ^ " exp 

physical parameters used in this work 



would be well protected by their planetary dynamos (if 
they actually exist or existed in the past) but that GJ 
667Cc has been too exposed to the erosive action of the 
stellar wind and probably has lost its atmosphere or its 
content of critical volatiles. 

In our conservative model we neglected the effect that 
an enhanced early stellar magnetic activity would have 
in the early erosion of the magnetically protected atmo- 
spheres, especially in the case of close-in planets around 
dM stars. It is interesting to consider how our results 
could be affected if a planet were exposed to more stress- 
ful conditions, e.g. coronal mass ejections and/or higher 
levels of XUV irradiation. We have performed sensitivity 
tests to study the effects of an enhanced stellar wind and 
larger fluxes of XUV radiation such as that expected for 
very active stars. To achieve this we modify the stellar 
wind pressure by maintaining the nominal velocity of the 
plasma (as predicted by the Parker's model) but increas- 
ing by a factor of 2 the density of stellar wind plasm£0. 
For the XUV flux we multiply by 2 the nominal XUV lu- 

^ Under typical conditions of solar CME, the velocity of the wind 
is not modified too much but the plasma densities scales up to 5-6 
times the average particle density 



minosity of G-K stars (Earth-mass planet) and maintain 
the initial high value of the XUV luminosity of dM stars 
(4.5 and 6.0 cases) simulating an extended period of 
stellar activity. The results of these sensitivity tests are 
depicted in the CME and Fxuv columns in figure [TT] 

CME conditions change the standoff radius is only 
10%. Again this could be explained by the weak depen- 
dence of this quantity on the stellar wind pressure and 
in particular to the number density of wind particles, eq. 

((SI) Rs ^ Psw^^ ~ n~^/^. A sustained increase in the 
XUV fluxes does not modify the standoff distance at all 
but it is able to noticeably increase the exobase radius 
and hence the time of exposure to the stellar wind ero- 
sion. As expected, the planet more susceptible to this 
effect is the Earth-mass planet which has a N/O-rich at- 
mosphere. Even under harsh conditions of stellar aggres- 
sion the mass-losses for TPs, irrespective of their mass, 
are no larger than 10 times the nominal value. Again this 
reinforces the conclusion that Gl 581d and HD 85512b 
are well protected by their potential PMFs. 

The existence of mobile lids is also a key factor for the 
existence of strong enough PMFs on habitable planets 
and it is required for the validity of the thermal evolution 
model results used in this work. The problem is being 
explored from different perspec tives and using com pli- 



mentary methodologies (s ee e.g. [Valencia et al.ll2007l and 
lO'NeiU fc Lenardidl2nn7l A consensus about the possi- 
bility of planets with greater mass having or not having 
lid activity has not yet been reached. However, different 
lines of evidence point to the fact that, in a wide range 
of rheological and thermal paramters, mobile lids could 
be c ommon in planets up to IOM0 (for recent results see 
e.g. iNoack et al.l[20Tl . 

To estimate the radius of the exobase our model re- 
lies on the results of detailed hydrodynamical and ther- 
modinamic models of highly irradiated atmospheres of 
TPs. However, we have used results calculated for two 
different types of atmosphere. Although in nature at- 
mospheric composition could also depend on planetary 
mass, for reliable comparisons it is better to use a sin- 
gle model. We also have a "void" in model results for 
planets in the mass range 1 — QMq. Although we have 
filled the 3.5 — 6 mass range using a simple linear 
extrapolation it is expected that future improvements to 
the atmospheric models will provide reliable values for 
the exosphere properties for planets in that interval of 
masses. We are confident however that the exobase ra- 
dius will not be too different from that used in this work 
and the main qualitative conclusions will not be modified 
significantly. 

6. SUMMARY AND CONCLUSIONS 

In the last few years we have seen significant improve- 
ments in the understanding of the role that planetary 
magnetic fields could have in the stability of the at- 
mopsheres of exoplanets. In line with these advances 
in this paper we have developed a comprehensive model 
of the evolution of magnetic protection of potentially 
habitable TPs, integrating in a single framework the 
results from very different specific areas of research in 
this field: thermal evolution of solid planets, scaling of 
dynamo-generated magnetic fields, magnetosphere mod- 
eling, physical properties of low-mass stars, stellar wind 
evolution and atmospheric modeling of highly irradiated 
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planets. 

Using this model we have studied the magnetic pro- 
tection of hypothetical habitable TPs in a wide range of 
planetary masses and have addressed for the first time 
the cases of the already discovered TPs found in the HZ 
of their host stars. In all these cases we have estimated 
the evolution of two key properties: the magnetopshere 
size as measured by the standoff distance and the ra- 
dius of the exobase. The direct comparison between 
these properties gives us information about the evolu- 
tion of the protection that the magnetosphere provides 
to the planetary atmosphere against the erosive action 
of the stellar wind. In order to estimate at to extent 
this magnetic protection prevents the loss of a signifi- 
cant fraction of the mass of the atmosphere or the loss of 
large amounts of critical volatiles (e.g. H2O and CO2) we 
have estimated the thermal-induced mass-loss. Knowing 
that non-thermal losses could be much larger, our conser- 
vative model provides an imderestimation of the stellar 
wind exposure effect. The potentially habitable TPs that 
under our model result in unsuitable conditions for mag- 
netic protection in reality are even worse than the model 
predicted. 

Our model is sensitive to a number of factors that af- 
fect the quantitative results we have obtained here. We 
have studied the sensitivity of our results to the expected 
variations in the model parameters and observed that the 
global conclusions are still very robust. 

We confirm that Earth-like planets, irrespective the 
composition of their atmospheres and even under the 
highest attainable dynamo-generated magnetic field 
strengths, will lose a significant fraction of their atmo- 
spheres or their critical volatile content if they are tidally 
locked in the HZ of dM stars. The case for the absence 
of habitable Earth-like planets around this type of star 
are almost closed. The case of habitable EPs around GK 
stars is not as good as previously expected either. Earth- 
mass planets with N/O-rich atmospheres, even under the 
best conditions of magnetic protection, will probably lose 
their atmospheres or their content of water if they are in 
HZ closer than ^ 0.8 AU. This limit excludes a large 
range of stellar masses (0.6-0.9 Mq) depending on the 
particular region inside the HZ (close to the inner or 
outer limit) where the planet resides. 

super-Earths with Mp > 3Mq seem to have better 
chances of preserving their atmospheres even if they are 
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tidally locked around dM stars. Under similar conditions 
of thermal and magnetic evolution there seems to exist a 
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